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INTRODUCTION 


If  a set  of  position -related  observations  are  made  simultaneously  on  a target  of 
interest,  the  joint  position  density  posterior  to  the  data  resulting  from  the  observations 
provides  all  of  the  information  available  and  the  basis  for  any  action  to  be  taken.  If  the 
observations  are  not  simultaneous,  or  if  some  time  has  elapsed  since  the  observations 
were  made,  and,  in  addition,  the  target  is  mobile,  the  introduction  of  a model  to  account 
for  target  motion  becomes  necessary. 

For  many  applications,  the  motion  may  be  conveniently  characterized  by  the  joint 
density  of  target  course  and  speed,  or,  more  precisely,  by  the  joint  density  of  target 
course  and  speed  during  the  particular  time  interval  of  interest.  If  the  characterization 
that  is  appropriate  for  this  time  interval  is  independent  of  time,  the  problem  of  specify- 
ing the  position  density  at  the  end  of  the  interval  is  greatly  simplified.  A treatment  in 
terms  of  space  variables  alone  is  then  possible  and  is  achieved  by  replacing  the  given 
joint  density  of  course  and  speed  by  that  of  course  and  range  traversed  in  the  time 
interval  that  is  applicable.  This  is  the  approach  used  in  reference  1 where  the  assump- 
tion of  stationarity  is  justified  by  the  application. 

In  many  cases,  however,  the  speed  or  course,  or  both,  will  change  during  the 
interval  of  interest.  The  interval  may  be  very  long,  the  target  may  be  confined  to  a 
limited  area  or  engaged  in  a maneuvering  tactic.  In  such  cases,  the  model  for  the 
motion  must  take  these  changes  into  account.  If,  however,  the  principal  concern  is 
the  final  position  of  the  target,  only  the  cumulative  effect  of  the  changes  which  take 
place  during  the  time  interval  contributes  to  this.  Thus,  if  the  vector  sum  of  all  the 
component  motions  can  be  determined,  the  problem  is  essentially  reduced  to  the  sta- 
tionary case  referred  to  above. 

For  the  application  at  present  under  consideration,  certain  assumptions  may  be 
made.  The  time  intervals  and  target  speeds  are  such  that  a flat -earth  approximation 
is  adequate.  The  components  of  the  motion  are  assumed  to  be  mutually  independent 
and  sufficiently  large  in  number  during  any  time  interval  that  will  occur  to  justify  the 
application  of  the  Central  Limit  Theorem.  Changes  in  course  will,  for  the  greater 
part  of  the  research  contribution,  be  uniformly  distributed.  The  principal  problem 
to  be  addressed  consists  of  finding  the  final  position  density  after  a time  interval  of 
arbitrary  length,  during  which  the  target  motion  is  of  a time -varying  nature,  given 
that  the  initial  position  is  normally  distributed.  In  addition  to  this  problem,  the 
research  contribution  provides  numerous  related  results.  All  of  these  results  are 
based  on  contributions  to  the  theory  of  wave  propagation  by  Rayleigh,  Rice,  Nakagami, 
and  others. 
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PRELIMINARY  REMARKS 

The  initial  position  of  the  target  is  designated  (x,  y)  and  the  joint  density  function 

p (x,  y)  is  assumed  known. 
rxy  } 


The  motion  of  the  target  during  any  time  interval  T is  assumed  to  consist  of  n 
segments  where  n is  large . An  arbitrary  segment,  say  the  jth,  is  characterized  by 


speed  V,  , course  cp  , and  duration  t.  or,  alternatively,  by  distance  traversed  or 

J J 


j 


In  addition  to  the  independence  of  different  segments,  the 


range  r.  = V.t.  and  course  cp. 

J J J J 

variables  and  r.  will  be  assumed  independent  of  cp. , and  the  latter  will  be  assumed  to  be 


uniformly  distributed  on  the  interval  0 cpj<  2tt  , for  all  j , throughout  most  of  the 
research  contribution.  Only  such  limitations  as  are  required  for  the  application  of 
tne  Central  Limit  Theorem  are  imposed  on  the  densities  p (r  ) and  p (V  ) . The 

rj  J j J 

projections  of  the  component  traverses  in  the  directions  of  the  reference  coordinate 
axes  are  given  by: 


x.  = r.  cos  cp,  = V.t,  cos  cp. 

J J j J j j 

y,  = r,  sin  cp.  = V.t.  sin  cp, 

j j '3  3 J j 


The  cartesian  and  polar  coordinates  at  the  end  of  the  time  interval  T are  denoted 
by  (u,v)  and  (R,  9)  respectively,  and  it  is  the  density  functions  of  these  variables  for 
which  formal  expressions  are  required.  The  relationship  of  these  variables  to  those 
previously  defined  is  illustrated  In  figure  1 . 


Note  that  (R,  0)  does  not  describe  the  closing  vector  of  the  random  motion  unless 
the  initial  position  coincides  with  the  origin  of  the  coordinate  system.  It  is  clear,  how 
ever,  that  if  an  additional  vector,  equivalent  in  effect  to  the  initial  conditions,  is  in- 
cluded as  part  of  the  motion,  the  final  position  density  will  be  unchanged,  but  (R,  0) 
will  then  describe  the  closing  vector  of  the  augmented  but  less  homogeneous  model  of 
the  motion  so  obtained. 


The  approach  that  will  be  taken  is  to  consider  a series  of  problems  of  increasing 
complexity  and  to  derive  expressions  relating  to  the  final  position  probabilities  in  each 
case.  These  expressions  will  contain  parameters  whose  values  may  be  ascertained 
from  those  used  to  characterize  the  random  motion.  For  the  application  of  interest, 
the  total  duration  of  the  motion  T will  be  known,  and  for  this  case  we  shall  assume  that 
the  mean-square  speed  V^3  is  also  specified.  This  is  defined  as: 

vTi  = i;0TEtva<t))dt  . 
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FIG.  1:  NOTATION 
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Since,  however, 


V(t)  = V.  for  t 

3 : 


<t  St 


j 


n 


this  may  also  be  written  as: 


V-E[f  EVTj] 


where  E denotes  expectation. 

We  shall  also  consider  the  case  in  which  T is  not  specified,  but  the  total  number 
of  segments  of  the  motion,  N , is  known.  While  this  case  has  no  immediate  bearing 
on  the  application  of  current  interest,  the  results  are  similar  in  form  involving  only  a 
redefinition  of  key  parameters.  For  this  case,  in  addition  to  N , we  shall  assume  that 
the  mean  square  range: 


is  specified. 

PROBLEM  1 

Problem  1 is  a special  case  of  the  problem  of  primary  interest.  The  initial  position 
is  known  and  coincides,  or  may  be  made  to  coincide,  with  the  origin  of  the  coordinate 


The  total  duration  of  the  motion  T and  the  mean  square  speed  are  specified. 


0 <r  cp^  < 2tt  . 
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n n 

U=E  Vi  COS(p  =£  X 

j=l  J J J j=l  J 


n 

= £ 
J=1 


n 


V)slntp)*j?1y) 


It  follows  from  these  relationships  that  the  variables  and  y.  are  also  identically 

distributed,  and,  since  n is  large,  the  conditions  of  the  Central  Limit  Theorem  are 
satisfied.  The  coordinates  of  the  final  position  u and  v are  therefore  approximately 
normal,  and  the  parameters  of  the  density  are  obtained  as  follows. 


For  any  time  t , 


9(0  = Sq  cp(t*)cit*  = £ cp  + 90 


J-l 


where  n'  is  the  integral  number  of  segments  traversed  in  the  time  t . But  since  the 
sum,  modulo  2tt  , of  any  number  of  random  variables  uniformly  distributed  on  the 
interval  0 £ cp^  < 2tt  is  uniformly  distributed  on  the  same  interval,  we  have 

0 i <p<t)  < 2tt  . 

Also,  since  V(t)  and  cp(t)  are  independent,  and 


u = V(t)  cos  cp(t)  dt 

v = /^V(t)  sincp(t)  dt 

E(u)=/0T  E [V(t)  ] E {cos  cp(t) } dt 
= 0 


since 


E {cos  cp(t) } = E {cos  cpj ) = ^ Jq2"1  cos  cp^dcp^  = 0 
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Similarly, 

E(v)  = 0 . 

au3  = E(u»)  = J*0T  J0TE{V(t)  V (f)}  E{cos  cp(t)  cos  <p(t’»dt  df 
= Tj'0TE{Va(t)}E{cosa9(t)}dt 
= I-/0T  E{Va(t)}dt 

since 

E{cosacp(t)}  = E{cosacp^}  =2^  JoTTcos8cpjdcpj  =?  * 
Therefore, 

°u’  = E<u*>  = T V - 

and,  similarly,  it  may  be  shown  that: 

cV*  - Bfr*)  - t V 

and  that: 

E(uv)  = 0 . 

Introducing  the  additional  notation: 

CT3  = — V 3 
CT  2 T 

<*  - 2aa  . 

we  obtain 

Pu(u)  = (2n)”^  a exp(-ua/2cra) 

Pv(v)  = (2n)"*  a exp(-va/2oa) 

Puv(u,v)  = pu(u)pv(v)  , 
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-i  -a 

(2tt)  ct  exp{ 


i (ua  +va)} 


Pr0(R,  0)  = R puy(R  cos  9,  R sin  0) 


0 < R < ® 


0 0 <:  2tt 


pe(e> 


(no)  J*  R exp(-Ra/o)dR 


-1 


(2tt) 


0 ^ 9 < 2rr 


-l 


2(a)  R exp  (-Ra/o) 


0 < R < « 


The  cartesian  coordinates  of  the  final  position  have  therefore  a bivariate  normal 
density.  The  marginal  density  of  0 is  uniform,  as  might  have  been  expected.  The 
distance  from  the  origin  to  the  final  position  is  Rayleigh  distributed  with  mean  square 
value  a . 

The  distribution  function  for  R is: 


pR(R)=  2(a)  JQR  u exp  (-u 3 /a)  du 
= 1 - exp(-Ra/or)  . 


Hence, 


Prob  (R  > p)  = exp(-pa/o) 


But,  since  a = E(Ra) 


and 
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Contours  of  constant  probability  are  concentric  circles  centered  at  the  origin.  In 
fact,  all  of  the  expressions  obtained  exhibit  this  central  symmetry  as  a consequence  of 
the  fact  that,  in  this  problem,  both  the  motion  and  the  initial  conditions  have  this 
property. 

PROBLEM  la 


This  problem  differs  from  Problem  1 only  in  the  following  respect.  Instead  of  the 
duration  of  the  motion  and  the  mean  square  speed,  the  total  number  of  segments,  N , 
and  the  mean  square  range  r^3  , are  specified. 


N N 

u = £ ri cos  % = £ x< 

j=i  J J j=i  J 


N 


v=£ 

j=i 


N 

r sin  cp  = £ y . 
J J j=i  J 


As  before,  the  conditions  of  the  Central  Limit  Theorem  are  satisfied.  Taking  expecta- 
tions, therefore,  we  obtain: 


E(u)  = E(v)  = 0 


E(ua)  = E(va)  = |e  (r3)  = ^ rNa 


E(uv)  = 0 . 

If,  therefore,  we  redefine  the  parameter  ga  for  this  problem  by: 


all  of  the  results  obtained  in  Problem  1 hold  for  this  problem  also. 

Further,  although  a characterization  of  the  motion  in  terms  of  its  duration  and 
mean  square  speed  will  be  assumed  for  the  remainder  of  the  problems  we  shall  consider, 
it  is  clear  from  the  above  that  formally  identical  solutions  will  hold  for  a characteriza- 
tion in  terms  of  total  number  of  segments  and  mean  square  range.  The  inclusion  of  this 
problem  is  to  indicate  the  single  modification  that  is  necessary. 
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PROBLEM  2 


The  statement  of  Problem  2 is  similar  to  that  of  Problem  1 with  one  exception. 

The  variables  r^  are  not  identically  distributed  for  all  j . 

In  the  treatment  of  the  previous  problems  it  was  possible  to  ascertain  the  densities 
of  final  position  variables  through  the  application  of  the  Central  Limit  Theorem  without 
reference  to  the  properties  of  the  density  p^  (r^)  . This  theorem,  however,  holds 

under  more  general  conditions  than  those  imposed,  so  all  of  the  results  obtained  re- 
main valid  when  the  variables  r,  have  different  distributions  for  different  values  of 

j , yet  satisfy  this  more  general  constraint. 


While  a precise  statement  of  the  most  general  conditions  under  which  the  Central 
Limit  Theorem  holds  is  available  (reference  2),  it  is  doubtful  if  it  would  serve  in  the 
present  context.  Essentially  what  is  required  is  that  the  variance  of  any  one  term  of 
the  process  is  negligible  in  comparison  to  the  variance  of  their  sum.  Thus,  if  the 


variables 
valid  if: 


r^  are  not  identically  distributed,  the  results  obtained  in  Problem  1 remain 


E(r^3)  = E(Vj8t.3)  «TaVTa 

for  all  j , and  in  Problem  la  if: 
E(r^)  « N rN3 

for  all  j . 


These  conditions  may  therefore  replace  the  requirement  that  the  component  ranges 
have  identical  distributions  in  Problems  1 and  la  respectively. 


To  complete  the  discussion  of  Problem  2,  we  consider  how  to  proceed  when  either 
of  these  conditions  is  not  satisfied.  If  the  condition  is  not  met  for  a single  value  of  j , 
the  contribution  of  this  component  may  be  separated  from  the  others  whose  combined 
effect  can  be  ascertained  via  the  Central  Limit  Theorem.  This  is  clearly  equivalent 
to  the  addition  of  a random  term  to  the  motion  described  in  Problem  1 . Motions  with 
this  feature  will  be  considered  later.  Where  the  condition  is  not  met  for  more  than  a 
single  value  of  j , the  sum  of  such  terms  may  be  treated  separately,  so  this  situation 
is  not  essentially  different. 
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PROBLEM  3 


In  Problem  3 the  initial  position  of  the  target  is  arbitrary  but  is  known  with  cer 
tainty.  All  other  factors  are  as  in  Problem  1 or  Problem  2. 

For  an  initial  position  (x,  y) 

rp 

u = x + J V(t)  cos  cp(t)  dt 

'o 

T 

v = y+  J V(t)  sin  cp(t)  dt  . 

0 

It  is  clear  from  the  preceding  problems  that: 

E(u)  = x 

E(v)  = y 

2 — 2—2 
au  - 'V  “ CT 

^ = . 

Therefore, 

Puv(u«v)  = Pu(u)pv(v) 

pRQR,  = (2tt)  a exp{-(u-x)2/2rr3  - (v-y)2/2a3} 

pRQ(R,  9)  = RPuy(R  cos  9,  R sin  0) 

= (it c)  R exp  £-  ^ fRs+r3-2Rr  cos  (9  -cp)  }J 

where 

r = (x3+y3)^ 
cp  = arctan  y/x  . 

pR(R)  = o_1  2R  exp{- ^(R3+ra)}  (2tt)  J*  ^expr^  Rr  cos  (0-cp) ) d0  , 
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= » -1  2R  exp  {-  — (Ra+ra) } (2n)_1  f ^ ‘ exp  f—  Rr  cos  0 ] d0 

-cp 

= o"1  2R  exp  { " ~ (R3+ra)  } 0<R<«  . 

This  is  the  Rice -Nakagami  Density. 

The  marginal  density  of  0 is  not  uniform  and  is  given  by: 


pe(0)  = r°pRe(R>e)dR 


= (2tt)  exp(-r3/a)  ^1  + Grr^ exp (Ga)  {l+erfG}J  , O50<2tt 


where 


G = a ^ r cos  (0  - cp)  . 

The  complement  of  the  distribution  function  of  R is  given  by: 

Prob  (R  > p)  = J*  p (u)  du 
P R 

= 2a  J " u exp  {-  ^ (u3+ra)  ) IQ  ( “■)  du  . 

Tables  of  this  integral  may  be  found  in  reference  3. 

The  solutions  to  this  problem  illustrate  the  asymmetry  introduced  by  a noncentral 
initial  position.  As  would  be  expected,  they  approach  the  corresponding  solutions  of 
Problems  1 or  la  as: 

=o  . 

RMS  o® 

PROBLEM  4 

In  this  problem  the  initial  position  is  uncertain  but  the  density  p (x,  y)  is  speci- 

xy 

fied.  Otherwise  the  motion  is  as  described  in  Problems  1 or  2. 
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From  Problem  3, 


Puv(u» v |xy>  = (2tt)  a"a  exp[_2^a  tu"x^a  + • 

Therefore, 

Puv(U,v)  = •TJ*  Puv(u,VlX’ y)  Pxy(x,y)  **  dy 

= (2tt)  _1a"aJ'J exp[-  ^ {(u-x)a  + (v-y)a}]pxy(x, y)  dx  dy  . 

If  (r,  cp)  are  the  polar  coordinates  of  the  initial  position  (x,  y)  , 
prcp(r»  9)  = rPxy(r  cos  «P»  r sin  9) 

Pr(r)  = ^o^Prcp^’  ^ dtP  * 

These  densities  may  therefore  be  determined  from  the  given  density  PXy(x»  y)  • 
Also  from  Problem  3, 

p^RlD-20-1  Rexp{-i  »•«■)]  . 

Therefore, 

Pr(r)  = Pr(r  |r)  Pr(r)  dr 

= 2a‘1  R exp  f-  i (Ra+ra) } Iq  (^)Pr<r)  dr 

and 

Prob(R  > p)  = 2*-1  J1"  J^“u  exp{-  | (u8+ra)}  ^ (~L)  pr(r)  dr  du  . 
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These  integral  expressions  permit  the  final  position  densities  v)  and 

p (R)  and  the  above  complement  of  the  final  position  distribution  P (R)  to  be  deter- 
K K 

mined  for  any  given  initial  position  density  p (x,  y)  , and  in  the  following  problem  we 

xy 

shall  obtain  explicit  results  for  the  case  where  the  initial  position  density  is  bivariate 
normal.  Before  considering  this  important  case,  however,  we  give  an  asymptotic  re- 
sult which  holds  for  any  initial  position  density. 

It  is  clear  that  for  sufficiently  large  values  of  R , the  influence  of  the  motion  is 
minimal  and  the  Prob(R>p)  is  largely  determined  by  the  initial  density  pr(r)  . In 

contrast,  for  small  values  of  R , the  Prob(R>p)  will  be  low  unless  r is  also  small. 
This  is  evident  from  the  exponential  term  in  the  integrand  of  the  above  expression  for 
this  probability.  But  if  r is  small,  the  asymmetry  is  small  and  PR(R)  may  be  approxi 
mated  by  the  Rayleigh  Distribution. 

A quantitative  treatment  of  these  issues  yields  the  following: 

pR(R)  = Pr(R)  for  R » afi 

= 2a"1  R exp(-  K3/a)  for  R «a^ 

Thus, 

Prob  (R  > p)  = j* " p (u)  du 
P r 

= 1 -Pr(p) 

= exp(-  p 3/a) 

PROBLEM  5 

Problem  5 is  the  special  case  of  Problem  4 where  the  initial  position  density 
Pxy(x,  y)  is  normal  with  parameters  px>  p^,  a3,  oy  3 and  p^  . 

Clearly  the  final  position  density  p (u,v)  may  be  obtained  from  the  relationship 
given  by  the  previous  problem  by  performing  the  necessary  integration.  It  is  obvious, 
however,  that,  for  this  case,  this  density  is  also  bivariate  normal,  and  that  its  param 
eters  p , p , cr  a,  c^8  and  puv  may  be  ascertained  directly  from  the  relationships: 


for  R»a“ 
for  R«a^ 
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U = X + u' 


v = y + v' 


v’  = E yi 

j J 


We  know  from  Problem  1 that  u'  and  v’  are  normal  variables  with  zero  means 
and  equal  variances  g3,  where  c3  may  be  obtained  from  the  parameters  which  charac- 
terize the  motion. 


Therefore, 


pu  = E(u)  = E(x)  = px 
Pv  = E(v)  = E(y)  = py  . 

Also,  since  the  motion  is  assumed  independent  of  the  initial  position, 

a _ 3,  3—  3.  3 

°u  ■ ax  + <V  - ax  + ° 


3 - _ 3 


V ■ ay*+  V*  ” ay3+  a 

P = a-1  a -1  cov («v) 

UV  U V 


= °u"1  av"1  [E(UV)  -E(u)E(v)] 

= au-1av-1  CE(xy)  - E(x)  E(y)] 

= ^ -1  a^1  cov  (xy) 

= au‘1av'1axOypXy 

_ CTxaypxy 

(axa  + as)i(ay3  + a3)£ 

This  completes  the  specification  of  the  final  position  density  p (u,  v)  , and  from  this, 
we  have: 


pRQ(R,  0)  = R Puy(R  cos  0,  R sin  0)  . 
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To  obtain  the  marginal  density  p (R)  it  is  convenient  to  rotate  the  coordinate 

R 

system  through  the  angle  3 where 


tan  23 


3 3" 

°u  % 


The  new  coordinates  are  then  given  in  terms  of  the  old  by  the  relations: 

u’  = u cos  0 + v sin  3 
vf  = -u  sin  3 + v cos  3 
and  the  old  in  terms  of  the  new  by: 

u = u’  cos  3 - v'  sin  3 
v = u'  sin  3 + v’  cos  3 . 

Also,  we  define  pu,  and  |iv,  by  the  relations: 

Li  , = u cos  3 - u sin  3 
*11  ru  v 

li  , = u sin  P + u cos  P 
'v  nu  v 

from  which, 

u = u , cos  p + u , sin  3 

Mi  U V 

li  = -li  , sin  P + u , cos  P . 

“v  ru  v 


Then 


where 


PuV.(u\v’)  = (2TT)^(Su,Sv,)'*exp  [- j {(u,-pu,)a/Su, 

•* £ ■■ fcr + v ’ I + ' 

\ 

Pr9(R,  9)  = R puV  {R  cos  (9-0),  R sin  (0-3))  . 


+ (V-Mv.)a/Sv,}] 
4(1~puv3>  |*1 

SiV  t ‘ ' 
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Thus, 

Q 

pr(r)  = jj)  R Pu,v,  fR  cos(fi-0),  R sin  (9-0)1  d0 

_ ^2^  0r  p fR  cos  Rsin0}d0 

-0 

= (2tt)  (Su,  Syt)  ^R  J2TT'^  exp  T-  j {(R  cos  0-pu,)  /Su,  + (R  sin  0-pv,)  /Sv,  }]d0. 

-0 

The  integral  on  the  right  of  this  expression  may  be  reduced  to  an  infinite  series  of  Modi- 
fied Bessel  Functions  of  the  First  Kind.  Details  of  the  reduction  are  similar  to  those  used 
in  the  appendix  of  reference  1 . When  this  reduction  is  carried  out,  we  obtain 

00 

p <R)  = (S  ,S  f)'*Rexp(-q)X;  (-l)kekIk(f)I2k(w)cos2kY  , 0<R<» 

k=0 

where 

= 1 for  k = 0 


= 2 for  k > 0 


This  is  the  Nakagami  Density  (reference  4). 
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Finally, 


P9(e)  = CPR9(R>9)dR 


R Puy(R  cos  0,  R sin  0 ) dR  . 


A closed  form  expression  for  this  density  has  been  obtained  only  in  special  cases. 
PROBLEM  5a 

The  statement  of  Problem  5a  is  identical  to  that  of  Problem  5.  The  difference  lies 
in  the  treatment. 

It  is  clear  that  part  of  the  complexity  of  the  results  obtained  in  the  preceding  prob- 
lem arises  from  the  fact  that  a prespecified  coordinate  system  was  assumed.  This,  of 
course,  is  necessary  if  any  information  relating  to  the  final  position  of  the  target  is  to 
be  conveyed  to  someone  else.  It  is  relatively  simple,  however,  to  devise  programs 
which  will  translate  and  rotate  the  coordinate  system  originally  specified  to  one  better 
adapted  to  the  problem  at  hand  and,  on  completion,  reestablish  the  original  reference. 
If  a program  of  this  type  is  available,  it  is  advantageous  to  center  the  new  coordinates 
at  the  mean  of  the  initial  position  density  (p  , u ) and  to  rotate  them  through  the  angle 
0 where  X y 


2(7  Q p 

t.n(3  = ^2-X£  _ 

°x  ~°y 

I 

Then,  if  (x\  y’)  denotes  the  initial  position  with  respect  to  the  new  coordinate  system, 
Px.  ,(x\y’)  = (2 tr)  (Sx,S  f^expj-  j:(x7Sx.  +y’ /S  )| 


where  S^S^,  are  given  by: 


x y xy  L x y l x y ux  y ' J 

Thus,  with  respect  to  the  modified  reference  system,  the  density  of  the  initial  position 
is  normal,  but  requires  only  two  parameters,  S , and  S , for  its  complete  specifica- 
tion. X 7 
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Continuing  as  in  Problem  5,  but  with  primed  variables  to  denote  reference  to  the 
modified  coordinate  system,  we  obtain: 

u , = u , =0 
ru  rv 

s,'  + «* 

V*  - sy,  + c3 

°u'V=0  ' 


Thus, 


PuV,(ut,V)  = pu,(u')pv,(V) 


= (2TTau,nv,)-1  exp  {-u*3/2au,s  -V3/2v3L 
PR,0,(R\  0,)=  R’PuV.  (R*  cose*,  R’sine*) 

R’3cos2e*  R,3sin29' 


= (2TTau,CTvl)-1  R*  exp 


'u* 


2V 


2tt 

pr,(r,)  - / pRt0t(R*,e*)de 

,,  r2Tt f R,s 

L" 


/TT  . K 2 d » 2 

= (2ttct  a.r  R ' f exp  (1+cos  29)  -7^(1- cos  29)  |d0 

u v 0 L ^°u*  av’ 


•]* 


= (c^.Oy,)-1  R*  exp 


R*3  / 1 1 \ 

4 (77‘  + ^~r 

u v J 

-R,scos  2 6 (—7--  T^r)  de 

V4°u*  4av*  'J 

- <au.av.)-‘  R-  exp  ["^r  (jTT*  + ^)j  >o| TT  ' ^)| 


(2n)_1  J21"  exp 
0 


0 < R'  < ® 


Curves  for  this  expression  and  its  integral  are  given  in  reference  5. 
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Pe.<9'>  = Jo”pR.0.<R,'e'>dR' 


* (2na„,<r  ,rl  J ’ R'  expl- 

u V o ' 


R,acosa0'  R,asina0, 


2V 


2v' 


dR' 


/o  W1  |cosa0*  , sina0'|_1  p®  , , 

I 25^7  +2^7"l  J0  “ P*P(->-a)  du 


— K 

2n  Kacosa0,+sina0' 


0 <;  0'  < 2tt 


where 


K = o ,/o  , 
u v 


Returning  to  the  joint  density  of  the  final  position,  it  is  readily  shown  that  curves 
of  constant  probability  are  ellipses  centered  at  the  new  origin  and  with  axes  parallel 
to  the  modified  coordinate  axes , Their  equations  are  given  by: 


for  values  of  C > 0 . 

The  parameter  C may  be  related  to  a measure  of  confidence  that  the  target  lies 
within  the  ellipse  which  it  defines.  If  this  measure  is  denoted  by  P , then 

C = -2  log(l-P)  . 

This  relationship  thus  defines  the  ellipse  corresponding  to  a prespecified  value  of  P . 
Its  area  is  given  by: 


TTau,av,C  = -2ttou,ctv,  log(l-P) 
and  the  lengths  of  the  semi-axes  by: 
ctu«  a°d  cty,  Jc~  . 

These  quantities  are  not  changed  when  the  original  coordinate  system  is  restored. 
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PROBLEM  6 


In  this  problem  the  initial  position  is  at  the  origin  of  the  coordinate  system  but 
changes  in  course  are  nonuniform.  The  components  of  the  motion  are  sufficiently 
large  in  number  and  homogeneous  in  composition  to  justify  the  application  of  the 
Central  Limit  Theorem.  As  before,  each  component  of  the  motion  is  assumed  to  be 
independent,  but,  at  least  initially,  we  shall  not  require  V.  and  co.  to  be  independent 
for  a common  value  of  j . J J 

As  before, 

T 

u = * V(t)  cos  cp(t)  dt  , 

0 

but 


T 

E(u)  f J*  E {V(t)  ] E {cos  cp(t)  ) dt 
0 

since  V(t)  and  cp(t)  are  not  independent.  Even  if  these  variables  were  independent 
for  all  t , since  cp^  » and  hence  cp(t)  , are  nonuniform,  in  general 

E {cos  cp(t)  } f 0 . 

Similar  arguments  prevail  for  E(v)  . 

Thus,  where  the  changes  in  course  are  nonuniform  but  the  conditions  of  the  Central 
Limit  Theorem  are  met,  the  final  position  density  p^u,  v)  is  normal  but  the  mean  is 

not,  in  general,  at  the  origin,  having  instead  some  value  (p  , p ) depending  on  the  joint 

densities  p„  (V.,cp.)  . 

V.cp.  j Y] 

The  variances  cua  and  are  obtained  by  taking  appropriate  expectations  but 

are  no  longer  expressible  in  terms  of  a parameter  of  the  speed  distribution  alone. 
Finally,  for  the  conditions  of  this  problem,  the  covariance  of  u and  v does  not 
automatically  vanish. 

The  joint  density  of  the  final  position  p^u.v)  is  thus  bivariate  normal  of  the  most 

general  type  requiring  the  five  parameters  p^,  p^,  au  2>  a^a>  and  for  its  complete 

specification.  This  being  the  case,  however,  all  of  the  results  of  Problem  5 are  also 
applicable  to  this  problem. 
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PR0<r,  0)  s R Puv(R  cos  0,  R sin  6) 


is  readily  obtained  once  the  parameters  of  puv(u,v)  have  been  determined. 

p (R)  has  a Nakagami  Density  with  parameters  derived  from  those  of  p (u,  v)  : 
R uv 


PR0(R,  6)  dR 


= R p (R  cos  0,  R sin  0)  dR  , 

•'O  uv 

As  previously  stated,  a closed-form  expression  for  this  density  is  not  available  in  the 
most  general  case.  There  is,  however,  a case  of  some  practical  interest  in  which  a 
simpler  expression  results  for  the  density  pR(R)  and  the  integral  on  the  right  of  the 
above  expression  can  be  evaluated. 


If  V.  and  cp, 
J j 


are  independent  for  all  j , and  cp.  is  symmetrically  distributed 


about  its  mean  value,  clearly  cp(t)  also  has  this  property. 


Therefore, 

r-p 

u =E(v)  = eJ  V(t)  sin  cp(t)  dt 
v 0 


= JTE{V(t)}E{sincp(t)}  dt 
0 


= 0 


since 


2tt 

E [sin  cp(t) ) = J sin  cp  p (cp)  dcp 

0 'P 


= J ^ sin  cp  p (cp)  dcp+  sin  cp  p (cp)  dcp 
0 * * 


= **  ^sincpp  (cp)dcp+  f°  sincpp  (cp)  dcp 
0 «P  p cp 

= 0 . 
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Therefore,  setting  My  = 0 in  the  Nakagami  Density,  we  obtain: 

pR(R)  = (SuSvf  * R exp(-q)  £ (-l)k  ^(OI^w)  , 0 < R < - 

k=l 


where 


fik  = 


1 for  k = 0 


= 2 for  k > 0 


«wg-+-(r+r) 

U U V 

™ 4(r-r) 

U V 


u 


w(R)  = R <r-  . 


u 


pe(9)  = tl+G'iexp(Gs)(l+erfG)]  , 0S8^2n 


where 


G(0)  = BK  cos  9 


1 + Ka 


2(K3cos30+sin30) 


B 

K 


u 


(a  + a ) 
x U V 

= a /a 
v u 
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